arXiv: 1505.00933v2 [astro-ph.HE] 19Jun2015 


Draft version June 22, 2015 

Preprint typeset using DT^X style emulateapj v. 5/2/11 


RECOLLIMATION SHOCKS IN MAGNETIZED RELATIVISTIC JETS 

Yosuke Mizuno 1 ' 2 , JoseL. Gomez 3 , Ken-Ichi Nishikawa 4 , Athina Meli 5 , Philip E. Hardee 6 , and Luciano Rezzolla 1 ' 7 

Draft version June 22, 2015 

ABSTRACT 

We have performed two-dimensional special-relativistic magnetohydrodynamic simulations of non¬ 
equilibrium over-pressured relativistic jets in cylindrical geometry. Multiple stationary recollimation shock 
and rarefaction structures are produced along the jet by the nonlinear interaction of shocks and rarefaction 
waves excited at the interface between the jet and the surrounding ambient medium. Although initially the jet 
is kinematically dominated, we have considered axial, toroidal and helical magnetic fields to investigate the 
effects of different magnetic-field topologies and strengths on the recollimation structures. We find that an ax¬ 
ial field introduces a larger effective gas-pressure and leads to stronger recollimation shocks and rarefactions, 
resulting in larger flow variations. The jet boost grows quadratically with the initial magnetic field. On the 
other hand, a toroidal field leads to weaker recollimation shocks and rarefactions, modifying significantly the 
jet structure after the first recollimation rarefaction and shock. The jet boost decreases systematically. For a 
helical field, instead, the behaviour depends on the magnetic pitch, with a phenomenology that ranges between 
the one seen for axial and toroidal magnetic fields, respectively. In general, however, a helical magnetic field 
yields a more complex shock and rarefaction substructure close to the inlet that significantly modifies the jet 
structure. The differences in shock structure resulting from different field configurations and strengths may 
have observable consequences for disturbances propagating through a stationary recollimation shock. 

Subject headings: galaxies: jets - shocks - magnetohydrodynamics (MHD) - methods: numerical 


1. INTRODUCTION 

Very Long Baseline Interferometry (VLBI) observations of 
Active Galactic Nuclei (AGN) jets often suggests the presence 
of quasi-stationary features (see, e.g., Jorstad et al. 2005; 
Lister et al. 2013; Cohen et al. 2014). These features can 
be associated with bends in the jet (see, e.g., Alberdi et al. 
2000 ) leading to enhanced emission due to differen¬ 
tial Doppler boosting (see, e.g., Gomez et al. 1993), 
or to recollimation shocks (see, e.g., Daly & Marscher 
1988; Gomez et al. 1995, 1997; Komissarov & Falle 1997; 
Casadio et al. 2013). Most of the observed quasi-stationary 
features appear in the innermost jet regions (Jorstad et al. 
2005, 2013; Fromm et al. 2013; Cohen et al. 2014) suggest¬ 
ing that they could be associated with recollimation, or recon¬ 
finement shocks produced by a pressure mismatch between 
the jet and the external medium. 

Multi-wavelength observations of Blazars suggest that 
high-energy 7 -ray flares are usually associated with the pass¬ 
ing of new superluminal components through the VLBI core, 
defined as the bright compact feature in the upstream end of 
the jet (see, e.g., Marscher et al. 2008, 2010). In order to pro¬ 
duce 7 -ray flares, an increase in particle and magnetic energy 
density is required when jet disturbances cross the radio core. 
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This increase in particle and magnetic density can naturally be 
explained by identifying the mm-VLBI radio core with a rec¬ 
ollimation shock (Gomez et al. 1995, 1997; Marscher et al. 
2010; Marscher 2014). Note however that at centimeter- 
wavelengths, the opacity core-shift observed in many sources 
(see, e.g., Kovalev et al. 2008; O’Sullivan & Gabuzda 2008; 
Sokolovsky et al. 2011) suggests that at these wavelengths 
the VLBI core corresponds to the transition between the opti¬ 
cally thick-thin regimes. 

Recollimation shocks have also been found at hundreds 
of parsecs from the central engine in 3C120, 3C380, 
and M 87 (Roca-Sogorb et al. 2010; Gabuzda et al. 2014; 
Asada & Nakamura 2012; Asada et al. 2014). The case of 
M87 is particularly interesting since the HST-1 complex, 
thought associated with a recollimation shock, shows behav¬ 
ior similar to a VLBI core, with new superluminal compo¬ 
nents emerging from its position (Giroletti et al. 2012). Very 
high-energy emission has also been observed in connection 
with variability in the HST-1 region (Cheung et al. 2007) 
similar to that observed in the VLBI core region of blazars. 

When a jet propagates through an ambient medium, pres¬ 
sure mismatch between the jet and the ambient medium 
naturally arises as a result of ambient pressure decrease. 
The pressure mismatch drives a radial oscillating motion 
of the jet and multiple recollimation regions inside the jet 
(see, e.g., Gomez et al. 1997; Komissarov & Lalle 1997; 
Agudo et al. 2001; Aloy et al. 2003; Roca-Sogorb et al. 
2008, 2009; Mimicaetal. 2009; Matsumoto et al. 2012). 
The resulting recollimation structure has been investi¬ 
gated analytically through a self-similar treatment of 
hydrodynamic jets (Kohler et al. 2012) and magnetized 
jets (Kohler & Begelman 2012), while Porth & Komissarov 
(2014) have studied the causality and stability of magnetized 
jets propagating in a decreasing pressure ambient medium. 
If a significant rarefaction wave is produced by the rec¬ 
ollimation and propagates into the jet interior, the ther- 
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mal energy of the plasma can be converted into kinetic en¬ 
ergy, increasing considerably the jet Lorentz factor. This 
is a purely relativistic effect, also referred to as the Aloy- 
Rezzolla (AR) booster (Aloy & Rezzolla 2006), which takes 
place in relativistic flows with a large tangential velocity 
discontinuity (Rezzolla & Zanotti 2002). Under these con¬ 
ditions, because the quantity 7 h is conserved across the 
rarefaction wave, with h the specific enthalpy and 7 the 
Lorentz factor, large jumps can take place in the latter, 
leading to a boost (Rezzolla & Zanotti 2013). This boost¬ 
ing mechanism is very basic and has been confirmed by a 
number of studies in hydrodynamical and magnetized jets 
(Mizuno et al. 2008; Komissarov et al. 2010; Zenitani et al. 
2010; Matsumoto et al. 2012; Sapountzis & Vlahakis 2013). 

In this paper we study in detail how various magnetic field 
configurations and strength affect the recollimation-shock 
structure using two-dimensional (2D) special relativistic mag¬ 
netohydrodynamic (RMHD) simulations. In particular, we fo¬ 


cus on the nonlinear rarefactions and shocks excited at the jet 
and ambient-medium boundary. Our present work provides 
an extension of the work presented by Gomez et al (1 997) to 
the case of a dynamically significant magnetic field. 

The paper is organized as follows. In Section 2 we describe 
the numerical method and setup used for our simulations. We 
present our results in Section 3, illustrating in detail the four 
different cases of purely hydrodynamical jets, as well as of 
jets endowed with axial, toroidal and helical magnetic fields. 
Finally we discuss the astrophysical implications in Section 
4, which also contains our conclusions. 

2. NUMERICAL SETUP 

We have perform 2D special-relativistic magnetohydrody¬ 
namics (SRMHD) simulations using the three-dimensional 
general relativistic MHD code RAISHIN adopting cylindrical 
coordinates (77, 6, z ) (Mizuno et al. 2006, 201 1). In particu¬ 
lar, we solve the SRMHD equations in the form 


9t ( 7 P) + ^d R (Rjpv R ) + -^<9, ( 7 pv*) + d z ( 7 pv z ) = 0, 
d t S R + ^ d R [7? (7 2 Hv r v r +p- b R b R )] + -^<9, (7 2 Hv r v* - b R b<+) + d z (7 2 Hv R v z - b R b z ) 


( 1 ) 

^Hv't’v 4 ’ +p- b 

R ’ 


( 2 ) 

dtS* + 4 d R [77 (7 2 Hv*v r - bh R )] + 4 < 9 , (7 2 77uV +p - b*b*) + d z ( 7 2 fTr;V - &V) = - b R b * , 

R R R 


d t S z + i d R [7? ( 1 2 Hv z v R - b z b R )] + (7 2 Hv z v* - b z b *) + d z (7 2 Hv z v z +p- b z b z ) = 0, (4) 

d t T + j^d R [77 (7 2 Hv r - b^ - 7 pv R )] + -^< 9 , (^Hv* - bH 1 - 7 -po*) d z (j 2 Hv z - b z R - 7/ cw z ) = 0, (5) 

dtB R + ^ ^ BR ~ ^B*) + d * ( ?*B R - vR B z ) = 0 . ( 6 ) 

+ -}-d r [7? {v r B^ - + d z ^B* - v^B z ) = vRb4> - v<t>BR J (7) 

d t B z + U [7? (v r B z - v z B r )1 + 4a, (v <t> B z - v z B^) = 0, (8) 

7? 77 


where we have set c = 1 and used Lorentz-Heaviside units. 
Here p is the rest-mass density, v is the plasma three-velocity, 
S is the three-momentum density, r is the conserved energy 
density, h = (e+p g )/p = l+e+p g /p is the specific enthalpy, 
with e = p{ 1 + e) the total energy density, e is the specific 
internal energy, and p g the gas pressure. Furthermore, B is 
the magnetic field measured in the Eulerian frame, while b a — 
(b*, b) is the magnetic field measured in the comoving frame, 
so that b 2 = b a b a = B 2 / 7 2 + {v ■ B) 2 , H = ph + b 2 is 
the total enthalpy, p m = b 2 / 2 is the magnetic pressure, and 
p = p g T Pm the total pressure. We also adopt as the equation 
of state that of an ideal gas with p g = pe{T — 1 ) and the 
adiabatic index V = 4/3 (Rezzolla & Zanotti 2013). 

In the RAISHIN code, a conservative high-resolution 
shock-capturing scheme is employed. The numerical fluxes 
are calculated using the Harten-Lax-van Leer (HLL) approx¬ 
imate Riemann solver (Hartenetal. 1983), and the flux- 
interpolated constrained transport (flux-CT) is used to main¬ 
tain a divergence-free magnetic field (Toth 2000). We use the 
monotonicity-preserving (MP5) spatial interoperation scheme 


(Suresh & Huynh 1997) and the third-order Runge-Kutta 
time-stepping scheme (Shu & Osher 1988) for all simula¬ 
tions. The RAISHIN code has the second-order accuracy 
based on flux-CT scheme even though we use higher spa¬ 
tial interoperation scheme (Mizuno et al. 2006). However, 
the higher-order spatial interoperation scheme leads to sharper 
transitions at the discontinuities. 

In the simulations presented here, a preexisting cylindri¬ 
cal flow is established across the simulation domain with 
jet radius Rj = 1. This setup represents a jet far be¬ 
hind the leading-edge Mach disk and bow shock (see, 
e.g., Mizuno et al. 2007, 2014). The flow is surrounded by 
a higher rest-mass density unmagnetized ambient medium. In 
all simulations the rest-mass density ratio is 77 = p g /p„ = 
5 x 10 :1 , where the subscripts j and a refer to jet and am¬ 
bient values, respectively. The ambient rest-mass density is 
constant with p a = 1.0 po- The jet speed is Vj = 0.9428 
in the z direction and 7 j = 3 with local Mach number 
M s = 1.69. We assume that the jet is initially uniformly over¬ 
pressured with p 3t j = 1.5 Pg^ a = l.bpgfi where p g $ is in 
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Figure 1. Case HD. 2D plots of: ( a ) the rest-mass density, ( b ) the gas pressure, and (c) the Lorentz factor. The white solid lines in panel {a) indicate the locations 
along the z axis used for the ID plots in Fig. 2. Shown instead in panel (c) with dashed, dotted and dash-dotted lines are the worldlines of the rarefaction waves 
(R), of the the shock waves ( S ), and of the contact discontinuities (C), respectively. 


units of po (see, e.g., Gomez et al. 1995, 1997; Agudo et al. 
2001; Mimica et al. 2009). The gas pressure in the ambient 
medium is constant. 

In order to investigate the effect of the magnetic field on 
the recollimation-shock structure, we have considered three 
different topologies: axial (or poloidal), toroidal, and helical. 
More specifically, for the axial-field case, the initial magnetic 
field B is uniform and parallel to the jet flow with B z = Bo, 
while for the toroidal-field case, we adopt the profile used by 
Lindetal. (1989) and Komissarov (1999) with 

f B Q R/R m if R < R m , 

B$ = \ B 0 R m /R if Rj> R >R m , (9) 

( 0 if R > Rj , 


where the magnetization radius R m = Rj /4. Finally, for 
the helical-field case, we chose a force-free helical mag¬ 
netic field with constant magnetic “pitch”, as the one used in 
Mizuno et al. (2009, 2011, 2014). The poloidal and toroidal 
magnetic field components are given by 


B z = 


B 0 


1 + ( R/a ) 2 


Bj, = 


Bo(R/a) 

1 + ( R/af 


( 10 ) 


where a is the characteristic radius of the magnetic field (the 
toroidal field component is maximum at radius a). The initial 
magnetic pitch is defined as 


= — (- — 
Rj Rj 


( 11 ) 


which is independent of Bo and such that a smaller Po refers 
to an increased magnetic helicity. In the magnetic helical-field 
case, we choose a = Rj / 2 , so that the initial magnetic field 
has constant helicity and pitch, with Po = 1/2. The values 
chosen for the initial Bq and the corresponding values of the 
plasma beta parameter, /3 p = p g /pm = 2 p g /b 2 (averaged 
and local minimum), as well as the magnetization parameter 
<t = b 2 /ph (averaged and local maximum) are listed in Table 


Case 

B-field 

Bo 

(dr) 

ftp, min 

(<r) 

0"max 

HD 

- 

0.0 

- 

- 

- 

- 

MHD-a 

axial 

0.1 

6.2 

5.3 

8.4 x 10~ 2 

9.0 x 10' 2 

MHD-b 

toroidal 

0.2 

140.5 

13.3 

1.1 x 10~ 2 

3.6 x 10 -2 

MHD-c 

helical 

0.2 

10.8 

1.3 

1.4 x 10' 1 

3.6 x IQ’ 1 


Table 1 

Basic properties of the various cases simulated. Listed are the magnetic-field 
topology, the initial magnetic-field strength Bo , the average and minimum 
plasma beta parameters (/3 P ) and /3p, m i n , as well as the average and 
maximum magnetization parameters (oj and cr max . 

1 for each model. Note that in all cases the relativistic jet is 
kinematically dominated for all values of Bq. 

The computational domain is 5 Rj x 30 Rj with a uniform 
grid of (R. z) = (128,300) computational zones. We impose 
outflow boundary conditions on the surfaces at R = J? max 
and z = 2 lri;ix - At z = 0 we use fixed boundary conditions 
and continuously inject the over-pressured jet into the compu¬ 
tational domain. The axisymmetry implies reflecting bound¬ 
ary conditions at R = 0. 


3. RESULTS 

In what follows we present the results of the simulations 
going through the four different magnetic-field configurations 
considered, i.e., cases HD, MHD-a, MHD-b, and MHD-c. 

3.1. Purely hydrodynamic jet 

We start our discussion of the results by illustrating what 
can be considered our reference configuration, that is, a purely 
hydrodynamicaljet (case HD). Figure 1 shows 2D plots of the 
rest-mass density, gas pressure and Lorentz factor for the hy¬ 
drodynamic case at t s = 200, where t s is in units of Rj. Mul¬ 
tiple recollimation shocks and rarefactions are evident along 
the jet propagation direction. Downstream of the inlet the 
over-pressured jet produces initially a weak conical shock that 
propagates into the ambient medium and a quasi-stationary 
conical rarefaction wave that propagates into the jet (dashed 
lines in panel (c) of Fig. 1). The initial shock is rather weak 
because of the small initial pressure discontinuity and leads to 
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Figure 2. Upper row: ID profiles perpendicular to jet axis of: ( a ) the rest-mass density, ( b ) the gas pressure, and (c) the Lorentz factor as measured at z = 4 
(red dotted), 7 (black solid), 14 (green dashed), and 20 (blue dash-dotted), respectively. Lower row: ID profiles along the jet axis (R = 0) of: ( d) the rest-mass 
density (solid) and the gas pressure (green dashed), and (e) the Lorentz factor. All panels refer to t s = 200. 


relatively smooth transition between jet and ambient medium. 
The over-pressured jet expands (following contact disconti¬ 
nuities, dash-dotted lines in Fig. 1) and the rarefaction wave 
leads to the conversion of thermal energy into kinetic energy 
of the jet (see, e.g., Aloy & Rezzolla 2006; Matsumoto et al. 
2012). The resulting acceleration is however rather modest 
given the choice made here for the initial conditions. 

One-dimensional (ID) profiles of hydrodynamic quantities 
perpendicular to the jet at axial positions z = 4, 7, 14, and 
20 are shown in the upper row of Fig. 2 and help interpret 
the dynamics of the recollimation shocks. Note that at z ~ 4 
the conical rarefaction wave converges to the jet axis and at 
this point the rest-mass density and pressure drop significantly 
and the flow is considerably accelerated, going from 7 ~ 3 
to ~ 4.6 (see the solid lines in Fig. 2). At the same time, 
the convergence of the rarefaction wave towards the jet axis 
leads to strong gas-pressure gradients, which slow jet expan¬ 
sion. The radial expansion of the jet ceases at z ~ 7 and at 
z>7 the jet starts to contract. Beyond z ~ 8 , a conical shock 
wave moves outward from the jet axis (see dotted lines in the 
panel (c) of Fig. 1) and reaches the jet edge at z ~ 13. At 
the conical shock, between z ~ 8 — 13, the rest-mass density 
and gas pressure increase and the jet Lorentz factor decreases. 
Furthermore, when the shock encounters the contact discon¬ 
tinuity at z ~ 13, the jet structure has mostly returned to the 
initial structure, albeit at slightly higher pressure and lower 
velocity in the jet. The innermost recollimation structures be¬ 
come stationary within several light crossing times. 

The panels in the bottom row of Fig. 2 offer a complemen¬ 
tary view and, in particular, they show ID profiles of hydro- 
dynamical quantities along the jet axis, i.e., at R = 0. Note 
that, as one would expect, the pressure profile follows exactly 
that of the rest-mass density as a result of the ideal-gas equa¬ 
tion of state employed. The Lorentz factor, on the other hand, 
is inversely proportional to the rest-mass density, so that the 
energy conversion from thermal to kinetic operated by the rar¬ 
efaction wave leads to a maximum in the Lorentz factor there 
where the rest-mass density has a minimum. Additional ac¬ 
celerations {i.e., peaks in the Lorentz factor) can be seen at 
z ~ 17 - 22, but these are less pronounced as the second 


set of rarefaction waves are not as efficient in reducing the 
rest-mass density locally. Note also that after the first cycle of 
recollimation shocks and rarefaction waves, the jet does not 
return to the initial conditions and so the subsequent recolli¬ 
mation shocks and rarefactions are weaker overall. 

3.2. Axial magnetic field 

Next, we consider the dynamics of the recollimation shocks 
when an axial magnetic field is present initially (case MHD-a). 
Figure 3 shows 2D plots of the rest-mass density, gas pressure, 
Lorentz factor and magnetic pressure for the axial field case 
with B 0 = 0.1 (MHD-a) at t s = 200. 

As in the hydrodynamic case discussed in the previous Sec¬ 
tion, immediately downstream of the inlet, the over-pressured 
axially magnetized jet produces a weak conical shock that 
propagates into the ambient medium and a quasi-stationary 
conical rarefaction wave that propagates into the jet. Also in 
this case, the flow is boosted as a result of the exchange from 
thermal energy to kinetic energy. More precisely, at z ~ 4.5 
the conical rarefaction wave converges to the jet axis and the 
rest-mass density, the gas- and the magnetic pressures drop 
significantly, leading to an acceleration of the flow from 7 ~ 3 
to ~ 6.0 (see the solid lines in Fig. 4). 

Note that the acceleration is larger than in the hydrody¬ 
namic case because the rarefaction wave is stronger and the 
latter is stronger because the axial magnetic field pressure acts 
in concert with the gas pressure but with different dependence 
on the jet radius. Similar results were obtained in 2D planar 
simulations of a flow bounded by a lower pressured ambi¬ 
ent medium (Mizuno et al. 2008; Zenitani et al. 2010). Fur¬ 
thermore, downstream of the rarefaction the gas pressure de¬ 
creases more than the magnetic pressure and the plasma beta 
decreases. 

Figure 4 reports ID profiles of various physical quantities 
perpendicular to the jet axis at axial positions z = 4, 7, 14, 
and 20 [see white solid lines in panel (a) of Fig. 3]. In this 
way it is possible to appreciate that the convergence of the 
rarefaction wave towards the jet axis leads to strong pressure 
gradients (both gas and magnetic), which slow jet expansion. 
Indeed, the radial expansion ceases at z ~ 7 and at z > 7 the 
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Figure 3. Case MHD-a. 2D plots of: (a) the rest-mass density, ( b ) the gas pressure, (c) the magnetic pressure, and (d) the Lorentz factor for the axial magnetic- 
field case with Bq = 0.1. The white solid lines in panel (a) indicate the locations along the z axis used for the ID plots in Fig. 4. All panels refer to t s = 200. 
The white solid lines in panel (a) indicate the ID plot locations in Fig. 4. 
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Figure 4. ID profiles perpendicular to jet axis of: (a) the rest-mass density, ( b ) the gas pressure, (c) the Lorentz factor, ( d) the axial magnetic field, ( e ) the 
magnetic pressure, and (f) the plasma beta. Different lines refer to 2 : = 4 (black solid), 7 (red dotted), 14 (green dashed), and 20 (blue dash-dotted). All 
quantities refer for the axial magnetic-field case (MHD-a) with Bq =0.1 and t s = 200. 


jet starts to contract. Beyond z cx. 8, a conical shock wave 
moves outward from the jet axis and reaches the jet edge at 
z ~ 13. When the shock encounters the contact discontinuity 
at z ~ 13, the jet has mostly returned to the initial structure, 
albeit at a slightly lower velocity. 

In summary, while maintaining many similarities with the 
purely hydrodynamic al evolution, the presence of an axial 
magnetic field leads to stronger recollimation shock and rar¬ 
efaction waves. In turn, because the AR booster is sensitive 
only to jumps in the specific enthalpy and not on whether a 
magnetic field is present, the sharper discontinuities in the 
flow lead to stronger boosts (see also the discussion in Sect. 

4). 

Finally, to assess the role played by the initial strength 
of the magnetic field Bq on the subsequent dynamics, we 
show in Fig. 5 the ID profiles of the gas pressure, of the 


Lorentz factor and of the magnetic-field strength along the jet 
axis, i.e., at R = 0. The different lines refer respectively to 
Bq = 0.0, i.e., the purely hydrodynamical evolution, as well 
as Bq = 0.025, 0.05, 0.075, and 0.1 respectively. The main 
feature to note in this case is that an increasingly stronger ax¬ 
ial magnetic field does not introduce qualitative changes in 
the flow dynamics. Indeed, the changes are only quantitative, 
with smaller values for the rest-mass density and pressure and 
consequently larger values of the Lorentz factor. It follows 
that the observational knowledge of the plasma velocity at dif¬ 
ferent positions in the jet could be used to deduce the strength 
of the magnetic field in the jet. 

3 . 3 . Toroidal magnetic field 

We continue our investigation by considering the dynamics 
of the recollimation shocks when a toroidal magnetic field is 
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Figure 5. ID profiles along the jet axis (R = 0 ) of: (a) the gas pressure, 
( b ) the Lorentz factor, and (c) the magnetic pressure for the axial magnetic- 
field case. Different lines refer to different values of the initial magnetic-field 
strength, i.e., Bo = 0 (black solid), 0.025 (red dotted), 0.05 (green dashed), 
0.075 (blue dash-dotted), and 0 . 1 , respectively (orange dash-double-dotted). 
All panels refer to t s = 200 . 

initially present (case MHD-b). In particular. Fig. 6 shows 2D 
plots of the rest-mass density, gas pressure, magnetic pressure 
and Lorentz factor for the toroidal field case with Bq = 0.2 at 
t s = 200. 

The general behavior is similar to both the hydrodynamic 
and axial magnetic-field cases. Namely, the over-pressured 
toroidally magnetized jet produces an initially weak conical 
shock that propagates into the ambient medium and, at the 
same time, a conical rarefaction wave that propagates into the 
jet. The latter expands radially, but the rarefaction waves pro¬ 
duced are also responsible for the conversion of the thermal 
energy of the plasma into kinetic energy of the jet following 
the basic mechanism of the AR booster. At z 4.5, this pro¬ 
cess leads to an increase in the Lorentz factor of flow that is 
accelerated from 7 ~ 3 to ~ 4.2 (see the solid lines in Fig. 
7c). 

Despite these analogies, the toroidal-field case also presents 
an important qualitative difference. Although it increases, the 
maximum value attained by the jet Lorentz factor is smaller 
than in the hydrodynamical case; this is in contrast with what 
was found for the axial magnetic-field case in the previous 
Section. The origin of this different behaviour is to be found 
in the magnetic tension introduced by the toroidal magnetic 
field and that essentially resists both inwards and outwards 
motions. This is quite visible in the behaviour along the 2 
direction of the magnetic pressure [see panel (c) of Fig. 6 ], 
which shows a “pinching” effect at z ~ 4,13 and z ~ 25, 
which also correspond to the locations where the rarefaction 
waves converge at the jet axis [see panel (a) of Fig. 6 ]. 

This structure in the magnetic pressure also appears in the 


radial profile of the plasma beta in Fig. 7, which reports ID 
profiles of several quantities perpendicular to the jet at axial 
positions z = 4, 7, 14, and 20. In particular, the panel (f) 
shows that plasma beta has two minima of the order of /3 p ~ 
20 on either side of the jet axis where the toroidal magnetic 
field is at a maximum (i.e., at II ~ ±0.3); this behaviour 
should be contrasted with the corresponding one shown in (f) 
of Fig. 4 for the axial magnetic-field case, where instead the 
smallest value of the beta parameter is fi p ~ 5. 

The presence of a toroidal magnetic field also alters con¬ 
siderably the flow after the shock encounters the contact dis¬ 
continuity at z ~ 13, breaking the appearance of a period¬ 
icity in the recollimation-shock structure. Interestingly, both 
the gas pressure and the Lorentz factor show an “0”-shaped 
structure around z ~ 16 — 23 [see panels (b) and (cl) of Fig. 

6 ] , with both of these quantities reaching a local minimum 
around z ~ 20 and 11 ~ 0 [see also panels (b) and (f) of Fig. 

7] . 

In summary, the presence of a toroidal magnetic field leads 
to weaker recollimation shocks and rarefaction waves, with a 
more complicated downstream structures than those found for 
the hydrodynamic and axial magnetic field cases. If observed, 
this different phenomenology should provide useful informa¬ 
tion to deduce the properties of the magnetic-field topology 
in the jet and, in particular, to establish whether this is purely 
toroidal. 

Finally, also in this case, we report in Fig. 8 the ID pro¬ 
files of the gas pressure and of the Lorentz factor along the 
jet axis, i.e., at II = 0 , as well as the magnetic-field strength 
along a direction slightly off the axis, i.e., at R = 0.24 (we 
recall that the magnetic field is zero along the axis). The fig¬ 
ure aims at establishing the variations introduced by the ini¬ 
tial magnetic-field strength Bq and so different lines refer to 
B 0 = 0.0 (i.e., the purely hydrodynamical case) as well as to 
Bq = 0.05, 0.1, 0.2, and 0.3. Perhaps not remarkably, only 
minimal changes are visible in the radial profiles of the rest- 
mass density, of the gas pressure and of the jet Lorentz factor 
as the magnetic field strength is increased. This is mostly due 
to the fact that the magnetic tension introduced by the toroidal 
magnetic field acts mostly in the R direction, producing only 
high-order variations in the direction of propagation of the jet. 
Such a behaviour is markedly different from the one encoun¬ 
tered in the case of an axial magnetic field, and indeed Fig. 8 
should be contrasted with the corresponding Fig. 5. 

3.4. Helical magnetic field 

Finally, we consider the dynamics of the recollimation 
shocks when a helical magnetic field is present initially (case 
MHD-c), with constant initial pitch P 0 = 1/2. As for the 
other cases, we show in Fig. 9 2D plots at t s = 200 of the 
rest-mass density, of the gas- and magnetic pressures, and of 
the Lorentz factor when Bo = 0.2. 

Similar to the other magnetic-field topologies, also in this 
case the over-pressured jet produces an initially weak conical 
shock that propagates into the ambient medium and a conical 
rarefaction wave that propagates into the jet. Again, conser¬ 
vation across the rarefaction wave of 7 h implies that a conver¬ 
sion of thermal energy to jet kinetic energy takes place across 
the rarefaction wave, with a consequent acceleration of the 
flow. Although the pitch considered here is less than one, so 
that the toroidal magnetic field is larger than the axial one at 
least initially, the effective behaviour of the plasma is closer to 
the one seen in the case of an axial magnetic field than in the 
case of a toroidal magnetic field. In particular, the increase in 
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Figure 6. Case MHD-b. The same as in Fig. 3 but for the toroidal magnetic-field case with Bq = 0.2. 
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Figure 7. The same as in Fig. 5, but for the toroidal-magnetic field case. 

the Lorentz factor is larger than in the hydrodynamical case, 
although only slightly, i.e., with downstream Lorentz factors 
going up to 7 ~ 5.0. On the other hand, as we have encoun¬ 
tered in the previous section for a purely toroidal magnetic 
field, the small pitch also implies that the plasma beta reach 
its minimum values there [cf, (f) in Fig. 7]. 

Figure 10 reports the ID profiles of several quantities per¬ 
pendicular to the jet axis at the positions z = 4, 7,14, and 20. 
Note the close analogies with the corresponding behaviours 
shown in Fig. 7 in the case of a purely toroidal magnetic field. 
Also here, the convergence of the rarefaction waves towards 
the jet axis leads to strong gas and magnetic-pressure gradi¬ 
ents, which ultimately slow the jet expansion. This ceases at 
z ~ 10 — 14, when weak shocks return the jet to an over¬ 
pressured structure and with a reduced Lorentz factor and 
rather different conditions than the initial ones. 

In summary, the presence of a helical magnetic field in the 
jet leads to a rather complex behaviour in both the recolli¬ 
mation shock and rarefaction structure. While the flow can 


be magnetically dominated with plasma beta of order unity, 
the Lorentz factor is increased with respect to the purely hy- 
drodynamical case, mimicking the behaviour seen for purely 
axial fields. This global behaviour is clearly influenced in 
part by our choice for the initial pitch and we expect that if a 
lower magnetic pitch parameter is chosen, corresponding to a 
more tightly wrapped helical field, the toroidal field will dom¬ 
inate over the axial field effects, with a sub-hydrodynamical 
Lorentz factor boost (see also the discussion in Section 4). 

Finally, the role played by the initial strength of the mag¬ 
netic field Bo on the subsequent dynamics is shown in Fig. 
11, which reports the ID profiles along the jet’s axis of the 
gas pressure, of the Lorentz factor and of the magnetic field 
strength. Again, the different lines refer respectively to Bq — 
0.0, and to B 0 = 0.05, 0.1, 0.15, and 0.2 (dotted, dashed, 
dash-dotted and dash-double-dotted lines, respectively). Also 
in this case, it is possible to observe how the phenomenol¬ 
ogy of the helical field is a combination of the one seen for 
purely axial and purely toroidal magnetic fields. In particu- 
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toroidal B, t=200 
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Figure 8. The same as in Fig. 5, but for the toroidal magnetic field case 
with Bo = 0 (black solid), 0.05 (red dotted), 0.1 (green dashed), 0.2 (blue 
dash-dotted), and 0.3 (orange dash-double-dotted). 

lar, as in the toroidal-field case, the gas pressure shows only 
modest changes along the z direction when the magnetic field 
is increased. On the other hand, as in the axial-field case, 
the Lorentz factor shows peaks that increase and move down¬ 
stream with increasing magnetic-field strengths. 

4. JET ACCELERATION 

In the previous Sections we have illustrated in detail how 
a robust feature of the recollimation shocks considered here 
is the acceleration of the flow downstream of the inlet as a 
result of the AR booster. We have also remarked that the oc¬ 
currence of this acceleration is the result of the conversion 
of the plasma thermal energy into kinetic energy of the jet, 
or, more precisely, of the conservation of the quantity 7 h 
across a rarefaction wave (Aloy & Rezzolla 2006). In addi¬ 
tion, we have discussed how different magnetic field topolo¬ 
gies lead to systematically different values of the maximum 
increase in the Lorentz factor measured downstream of the 
inlet, i.e., A 7 max = max( 7 ) — 70 . What we have not yet dis¬ 
cussed, however is how this accelerating mechanism depends 
on our choice of the initial magnetic-field strength, Bq . 


This point is addressed in the left panel of Fig. 12, 


which reports the relative increase of the Lorentz factor 
with respect to the purely hydrodynamic al evolution, namely, 
A7max/(7max)nD = 7max/(7max)nD 1, as a function of 
the initial magnetic field. Obviously, this quantity can either 
be positive or negative and provides a direct measure of the 
fractional boost. Shown with different symbols are the dif¬ 
ferent magnetic-field topologies, with crosses referring to the 
axial magnetic field (case MHD-a), diamonds to the helical 
magnetic-field (case MHD-c), and star crosses to the toroidal 
magnetic-field (case MHD-b). 

When presented in this manner, it is then straightforward 
to realize that axial and helical initial magnetic fields lead 
to Lorentz boosts that are larger than in the hydrodynami- 
cal case, while the opposite is true for purely toroidal mag¬ 
netic fields, for which an acceleration is still present but this 
is smaller than in the hydrodynamical case. We have already 
discussed in the previous Sections that the origin of this dif¬ 
ferent behaviour has to be found in the fact that an axial mag¬ 
netic field adds an effective gas pressure and results in larger 
rest-mass density and pressure gradients across the rarefac¬ 
tion waves induced downstream of the inlet. In turn, the AR 
booster translates these stronger waves into larger accelera¬ 
tions of the flow. It is also instructive that, in the case of 
a purely axial magnetic field, the behaviour of the relative 
boost has a simple quadratic dependence on the initial mag¬ 
netic field. This is simply because the relative boost scales 
as 

A7 ma x P - Pm j-,2 

7 - 1 —OC-1 = — oc B z . (12) 

(7max)HD Pg Pg 

This is confirmed by the very good match between the nu¬ 
merical data and a quadratic fit, which is indicated with a red 
dashed line. 

Also shown in Fig. 12, but in the middle panel, are the 
fractional differences in the specific enthalpy at the position 
of the maximum Lorentz factor, i.e., h/(h) un — 1. Here too, 
different symbols refer to the various magnetic-field topolo¬ 
gies ( cf., left panel) and different initial strengths. Clearly, 
this panel offers a complementary view to the one in the left 
panel and summarizes much of what already discussed in the 
previous Sections. Namely, that axial and helical magnetic 
fields leads to smaller values of the specific enthalpy in the 
downstream solution, in contrast to the case of toroidal mag¬ 
netic fields where instead the values of the specific enthalpy 
at Lorentz-factor maximum increases with the initial mag¬ 
netic field. Finally, shown in the right panel of Fig. 12 is 
the behaviour of the minimum plasma beta /3 p , m in as a func¬ 
tion of the initial magnetic field. Interestingly, all magnetic- 
field topologies show the same (and expected) quadratic de¬ 
pendence as Bq 2 . 

4.1. Dependence on magnetic pitch 

As discussed in Section 2, in the case of a helical magnetic 
field, we have an additional degree of freedom represented by 
the initial magnetic pitch as defined in Eq. (11). By suitably 
choosing the initial pitch, i.e., the ratio a/Rj, it is possible 
to scan the range of possible magnetic-field configurations, 
which range from an essentially toroidal magnetic field' for 
P 0 < 1 to an axial axial magnetic field for Pq 1 . 


8 Note that in practice even if the magnetic pitch is very small, the magnetic 
field does not reach the toroidal magnetic-field profile discussed in Section 
3.3 because there is always a nonzero poloidal component at the jet center 
R = 0. 












Recollimation Shocks in Magnetized Relativistic Jets 


9 


helical B (MHD-c), t=200, B o =0.2 
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Figure 9. Case MHD-c. The same as in Fig. 3, but for the helical magnetic-field case with Bq = 0.2. 


helical B (MHD-c), t=200 , B 0 =0.2 

(a) p (b) 



Pg 




R 


R 


(c) Y 



R 

(f) Pp 



R 


Figure 10. ID profiles perpendicular to the jet axis of: (a) the rest-mass density, (b ) the gas pressure, (c) the Lorentz factor, ( d) the toroidal magnetic field, ( e) 
the axial magnetic field, and (/) the plasma beta as measure at z = 4 (black solid), 7 (red dotted), 14 (green dashed), and 20 (blue dash-dotted). All panels refer 
to the helical magnetic-field case with B o = 0.2 at t s = 200. 


Given this freedom in the parameterization, and given that 
different magnetic-field topologies lead to different amplifi¬ 
cations to the AR booster, it is interesting to ask whether a 
correlation exists between the initial pitch and the increase in 
the maximum Lorentz factor downstream of the inlet. This is 
shown in Fig. 13, which shows A 7 max /( 7 max ) HD as a func¬ 
tion of the initial pitch l\, for a number of simulations carried 
out with Bq = 0.1. Interestingly, we found that the rela¬ 
tive increase in the maximum Lorentz factor has a very clear 
dependence with the pitch, smoothly joining the two extreme 
cases of a toroidal and axial magnetic fields, respectively. Fur¬ 
thermore, the transition between the two regimes takes place 
at Pq > 1, that is, when a > Rj, saturating to the axial field 
case when a ~ 10 Rj. Finally, the dependence can be accu¬ 
rately approximated with a simple expression of the type 

'Tmax — Cl + c 2 tanh[c 3 (Po - 1.0)], (13) 

where ci ~ 5.43, c 2 — 0.6, and c 3 ~ 0.8, and is indicated 


with the red dashed line in Fig. 13. 

The relevance of expression (13) is that it provides, at least 
in principle, yet another useful tool to deduce the properties 
of the magnetic-field topology from the observations of the jet 
dynamics. Indeed, if radio-astronomical observations could 
provide a reliable measure of the Lorentz factor in the bulk of 
the jet and at its base, i.e., 7 max and 7 j, then expression (13) 
would provide a simple way to deduce the degree of helicity 
of the magnetic field in the jet. This is because while the re¬ 
lation between 7 max and the pitch Pq changes with the initial 
Lorentz factor in the jet 7 j, the fact that 7 max and 7 j scale lin¬ 
early [cf, left panel of Fig. 15] implies that the functional de¬ 
pendence of 7 max / 7 j with Pq will not change. This is shown 
in the right panel of Fig. 13. Hence, once and if 7 max/ 7 j is 
measured, a direct estimate on P 0 is possible, which is only 
weakly dependent on the jet overpressure. Conversely, if the 
magnetic field pitch angle could be determined, e.g., from po- 
larimetric observations, then it could be possible to estimate 
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Figure 11. The same as in Fig. 5, but for the helical magnetic-field case 
with Bo = 0 (black solid), 0.05 (red dotted), 0.1 (green dashed), 0.15 (blue 
dash-dotted), and 0.2 (orange dash-double-dotted). 

7max- Such a measurement would obviously be of great im¬ 
portance for the interpretation of the very rapid TeV variabil¬ 
ity in AGN jets. 

5. CONCLUSION 

We have performed 2D SRMHD simulations of the prop¬ 
agation of an over-pressured relativistic jet leading to the 
generation of a series of recollimation shocks and rarefac¬ 
tion waves. The overall dynamics of this process is rather 
well known. Downstream of the inlet, the jet produces a 
weak conical shock that propagates into the ambient medium 
and a quasi-stationary conical rarefaction wave that propa¬ 
gates inwards. The significant drops in rest-mass density 
and pressure produced by the rarefaction waves are also re¬ 
sponsible for the conversion of the thermal energy into ki¬ 
netic energy of the jet (see, e.g., Aloy & Rezzolla 2006; 
Matsumoto et al. 2012). This is a purely relativistic ef¬ 
fect that develops in the presence of strong tangential dis¬ 
continuities in the flow (Rezzolla & Zanotti 2002), and that 
leads to a consistent and robust boost of the fluid. In ad¬ 
dition, the nonlinear interaction of the shocks and rarefac¬ 
tion waves lead to a stationary and multiple recollimation- 
shock structure along the jet, which has been reproduced 
by a number of the hydrodynamic simulations under dif¬ 
ferent physical conditions (see, e.g., Gomez et al. 1997; 
Komissarov & Falle 1997; Agudo et al. 2001; Aloy et al. 
2003; Roca-Sogorb et al. 2008, 2009; Mimica et al. 2009; 
Matsumoto et al. 2012). 

We have here extended previous hydrodynamic work to de¬ 
termine the effect that magnetic fields with axial, toroidal, 
and helical topologies have on the jet dynamics. In partic¬ 
ular, we have found that an axial magnetic field behaves as 
an additional gas pressure, leading to a sharper recollimation- 


shock structure and thus larger flow accelerations. On the 
other hand, a toroidal magnetic field tends to reduce the in¬ 
ward motion of the rarefaction waves, thus leading to a weaker 
recollimation-shock structure and accelerations that are even 
smaller than in the case of a purely hydrodynamical evolu¬ 
tion. Finally, a helical magnetic field tends to yield a be¬ 
haviour which is a combination of those observed for purely 
axial/toroidal magnetic fields, with larger Lorentz factors, but 
also with a rather complex recollimation-shock substructure. 

As predicted by the basic properties of the AR booster, 
we have found that the boost in the jet flow, as measured in 
terms of the maximum Lorentz factor, is anti-correlated with 
the values of the specific enthalpy. Furthermore, the boost 
grows quadratically with the initial magnetic field in the case 
of a purely axial flow, at smaller rates in the case of a helical 
magnetic field, and it decreases systematically when a purely 
toroidal magnetic field is present. Finally, we have shown that 
the maximum Lorentz factor exhibits a smooth behaviour in 
terms of the initial pitch in the case of an axial magnetic field. 

Stationary components are commonly seen in parsec-scale 
VLBI observations of AGN jets (see, e.g., Jorstad et al. 
2005; Lister et al. 2013; Cohen et al. 201 4), being normally 
associated with recollimation shocks. Studying in detail the 
structure in these stationary components for a direct compar¬ 
ison with our simulations requires however resolving the jet 
structure across the jet width. This can be achieved through 
“regular” cm-VLBI observations of nearby sources, like M87, 
but for the majority of the AGN jets it would require a sig¬ 
nificant increase in angular resolution. This can be obtained 
through VLBI observations at even shorter wavelengths, in 
which the 3 mm GMVA (Global Millimeter VLBI Array) and 
1.3 mm VLBI observations with the Event Horizon Telescope 
and Black-Hole Cam projects can achieve angular resolutions 
between 50 and 20 microarcseconds (see, e.g.. Fish et al. 
2014; Krichbaum et al. 2013), and through space VLBI ob¬ 
servations. The space VLBI mission “RadioAstron” has re¬ 
cently successfully achieved ground-space fringe detections 
for observations at 1.3 cm on baselines longer than eight Earth 
diameters, allowing imaging the innermost regions in AGN 
jets with an unprecedented angular resolution of ~ 20 mi¬ 
croarcseconds (Gomez et al. 2015). 

Whilst the range of simulations carried out here have tried 
to explore a rather large portion of the space of parame¬ 
ters, a number of improvements on the approach followed 
here can be made. First, even though we have assumed 
an over-pressured jet to produce the recollimation-shock 
structure, a pressure mismatch between jet and ambient 
medium arises naturally if the pressure in the ambient 
medium decreases with distance from the central object. In 
this case, the recollimation-shock structure depends on the 
ambient-medium pressure scale (see, e.g., Gomez et al. 
1997; Komissarov & Falle 1997; Agudo et al. 2001; 
Aloy et al. 2003; Mimica et al. 2009; Kohler et al. 2012; 
Matsumoto et al. 2012; Porth & Komissarov 2014). Al¬ 
though we do not expect that significant qualitative differ¬ 
ences will emerge from a more realistic modelling of the 
ambient medium, we will extend the current investigation 
to include larger overpressure ratios, as well as declining 
pressure profiles outside the jet. 

Second, in this work we assume the relativistic jet is kine¬ 
matically dominated initially for all values of the initial mag¬ 
netic field. However, magnetically dominated relativistic jets 
(Poynting-flux dominated jets) might be plausible from the re¬ 
sults of previous analytical (e.g., Heyvaerts & Norman 2003; 
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Figure 12. Dependence of: (a) the relative difference of the maximum Lorentz factor 7max/(7max) H D — L (O the relative difference of the specific enthalpy 
h/ (/i) IID — 1, and (c) the minimum plasma beta. All quantities are shown as a function of the initial magnetic-field strength. Bo in the case of axial (red crosses), 
toroidal (blue stars), and helical magnetic fields (black diamonds) with Po = 0.5. The red dashed line in panel ( a ) indicates a quadratic fit [cfi, Eq. (12)], while 
the black dot-dashed line in panel (c) refers to a Bq scaling. 





Figure 13. Left panel: Dependence of the maximum Lorentz factor 7 max (right vertical axis) and relative difference of the maximum Lorentz factor relative 
to the purely hydrodynamical case 7max/(7max) HD — 1 (left vertical axis) as a function of the initial magnetic pitch, Po; the initial jet Lorenz factor is set to 
7 j = 3; the red dashed line indicates the fitting with a hyperbolic tangent function [cf, Eq. (13)]. Central panel: Dependence of the maximum Lorentz factor 
as a function of the initial magnetic pitch, for three different values of the initial jet Lorenz factor 7 j = 2 (stars), 3 (diamonds), and 4 (triangles); the red dashed 
line indicates the fitting as in the left panel. Right panel:. The same as in the central panel but when the maximum Lorentz factors are normalized to the initial 
one; note that in this case all curves essentially overlap. All cases refer to an initial magnetic field Bq = 0.1. 


Beskin & Nokhrina 2009; Lyubarsky 2009) and numerical 
studies (e.g., Komissarov et al. 2006; Tchekhovskoy et al. 
2009; Porth et al. 2011; McKinney et al. 2012). We will in¬ 
vestigate this type of relativistic jets in future work. 

Third, it is clear that the magnetic modifications of the 
recollimation-shock structure found here may affect in a sig¬ 
nificant way the signatures of stationary components seen in 
several AGN jets, especially in polarized flux. An even clearer 
signature difference may be obtained when a shock like 
disturbance propagates through the stationary recollimation- 
shock structure. In view of this, we will investigate the propa¬ 
gation of such a shock-like disturbance through a magnetized 
recollimation shock and calculate the corresponding emis¬ 
sion following the approach presented in, e.g., Gomez et al. 
(1997). 

Finally, although our simulations show that the 
recollimation-shock strength depends on the jet magne¬ 
tization and magnetic field structure, it is clear that its 
observed appearance as imaged by VLBI would also depend 
on the jet-viewing angle, and is obviously influenced by 
Doppler-boosting effects. Hence, depending on the jet 
bulk flow Lorentz factor and viewing angle, the Doppler 
boosting may cancel, or even reverse, the increased emis- 
sivity obtained in the recollimation shock and due to the 
enhanced particle and magnetic-field energy rest-mass 


density (see, e.g., Gomez et al. 1997; Aloyetal. 2003; 
Roca-Sogorb et al. 2008, 2009; Mimica et al. 2009). These 
effects will also be the focus of future investigations. 
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HD, t=600, B 0 =0.0 
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Figure 14. Upper panels: 2D plots of: (a) the rest-mass density, (b) the gas pressure, and (c) the Lorentz factor. Lower panels: ID profiles along the jet 
axis (R = 0) of: (d) the rest-mass density (solid) and the gas pressure (green dashed), and (e) the Lorentz factor. All panels refer to t 3 = 600 for a purely 
hydrodynamical jet and an extent in the z direction that is three times larger than the one considered in the main text. 


APPENDIX 

CONVERGENCE TESTS 
Large-scale hydrodynamical jet 

Due to the uniform ambient medium, the modifications of the recollimation-shock and rarefaction-wave structures is expected 
to be rather small and a quasi-periodic structure should develop on larger axial scales (see also Matsumoto et al. 2012). In 
this Appendix we test this assumption, and hence the development of a self-similarity in the recollimation-shock structure, by 
considering a purely hydrodynamical simulation having the same resolution discussed in Section 3.1, but with an extent in the z 
direction from 0 to 90, that is, three times larger than what presented in the main text. The results of this simulation are reported 
in Fig. 14, with the top three panels referring to a two-dimensional view, and the bottom two panels to a cut along the z axis 
at t s = 600. It is clear that the shock/rarefaction structure is periodic and with periodicity ~ 12 Ii 3 . This separation distance 
depends on the initial jet velocity, over-pressure ratio between jet and ambient pressure, and opening angle of the jet. 
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Figure 15. Dependence of the maximum Lorentz factor in a purely hydrodynamical jet (7y e t) H D as a function of: (a) the initial jet Lorentz factor "[ with fixed 
over-pressure ratio Pg,j/Pg,a = 1.5 and (b) the over-pressure ratio Pg,j/Pg,a with fixed initial jet Lorentz factor 7j = 3. The red dashed lines indicate the 
linear fitting. 

Dependence on initial jet Lorentz factor and over-pressure ratio of the maximum jet Lorentz factor 

As discussed in Section 4, the relative difference of the maximum jet Lorentz factor with respect to the purely hydrodynamic 
case depends on the initial magnetic field strength and magnetic pitch. It is important to know the maximum jet Lorentz factor of 
a pure hydrodynamic case. In this Appendix, we have investigated the dependence of the maximum jet Lorentz factor of a pure 
hydrodynamic jet on the initial jet Lorentz factor and over-pressure ratio. Figure 15 shows f/rnax )un as a function of the initial 
jet Lorentz factor ( 7 ,-) and of the over-pressure ratio (Pg,j /Pg,a)- Clearly, (7 jet ) TTD has a simple linear dependence on both the 
initial jet Lorentz factor and the over-pressure ratio with a coefficient of ~ 1.8 and ~ 3.75 respectively, which we indicate with 
the red dashed lines in Fig. 15. 
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